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Abstract 


The  concept  of  a  locally  stationary  random  process  is 
introduced  and  the  harmonic  analyis  of  such  a  process  is  examined, 
the  generalized  Wiener-Khintchine  relations  obeyed  by  locally 
stationary  random  processes  are  discussed.  The  dielectric  tur- 
bulence responsible  for  beyond-the-horizon  scattering  of  radio 
waves  is  cited  as  an  example  of  a  locally  stationary  random 
process. 
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1.  Preliminary  considerations 

In  this  note  we  introduce  the  concept  of  a  locally  stationary  random 
process,  which  is  a  natural  generalization  of  the  familiar  notion  of  a  station- 
ary random  process.  We  shall  exhibit  a  remarkable  symmetry  between  the  co- 
variance  and  the  spectral  density  of  a  locally  stationary  random  process,  which 
gives  insight  into  the  structure  of  stationary  random  processes.  First  of  all, 
we  recall  some  relevant  facts  about  random  processes  in  general,  and  stationary 
random  processes  in  particular  *-.  J  »  L  J  # 

Let  x(t)  be  a  random  process  (in  general  complex)  with  mean  zero. 
Then  by  the  covariance  of  x(t)  is  meant  the  quantity 

(1)  P(t,f)  =  <x(t)x*(t')>   , 

where  the  angular  parentheses  denote  the  ensemble  average  and  the  asterisk  de- 
notes the  complex  conjugate.   It  is  apparent  from  (l)  that  J   (t,f)  =  |    (t',t), 
Suppose  x(t)  is  harmonizable,  i.e.,  that  x(t)  can  be  written  as  the  Fourier- 
Stieltjes  transform 

(2)  x(t)  =    exp(icot)dg(u) 


■/ 


of  another  random  process  g(w),  where  the  integral  is  to  be  taken  in  the  mean- 
square  sense.  (Here,  as  below,  all  integrals  are  taken  from  -oo  to  oo.)  Then 
it  follows  that  the  covariance  of  x(t)  can  be  written  as 


(3) 


P(t,t»)      =  exp[i(ayt  -  co't')]  T(co,co»)   dco  dw« 


-  2  - 


where  the  two-dimensional  spectral  density  Y(w,w')  is  related  to  the  increments 
of  the  random  process  g(u)  by 

(U)  Y(co,co»)dco  dco'   =  <dg(u)-dg*(o)')>   ♦ 

In  writing  (3)  and  (ij),  we  assume  the  physical  case  of  absolutely  continuous 
spectral  densities.  Moreover  we  have 


(S) 


T(co,co')  =  —i-g.       exp[-i(wt- co't')]  P(t,t')dtdt'  , 


so  that  I   (t,f)  and  1(00,00')  are  two-dimensional  Courier  transform  pairs. 

The  zero-mean  random  process  x(t)  is  said  to  be  stationary  (in  the 
wide  sense)  if  its  covariance  is  of  the  form 

P(t,t')   =  <|x|2>  R(t-t')  . 

The  quantity  <^|x|  ^  is  the  average  instantaneous  power  of  the  process  x(t) 
and  is  independent  of  time;  R(t-t')  is  the  normalized  correlation  function  of 
x(t),  with  R(0)  ■  1.  As  is  well-known,  in  the  absolutely  continuous  case 
R(t-t')  is  the  Fourier  transform  of  a  non-negative  function  (|  (00),  the  power 
spectral  density  of  x(t),  i.e. 


(6) 
(7) 


R(r)  =      exp(iooT)  J  (co)  dco  , 

I   (w)  =   ^_     exp(-icoT)  R(r)  &r    . 
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Equations  (6)  and  (7),  together  with  the  requirements  on  R(T)  and  <$  (co), 
namely  that  R(T)   be  a  stationary  covariance  or  correlation  function  and 
that  <|  (co)  be  non-negative,  are  the  celebrated  Wiener-Khintchine  relations. 
It  is  to  be  noted  that  <|  (co)  is  normalized,  i.e., 


/ 


^  (co)  doo  =  1  • 


The  harmonic  decomposition  theorem  for  stationary  random  processes 
consists  of  the  statement  that  the  relation  between  the  increments  of  the  ran- 
dom process  g(co)  appearing  in  (2)  and  <|  (co)  is 

<dg(co)dg*(co«)>   =  <|x|  >  jj)  (co)6(co-to')dco  dec'   , 
or  equivalently 

(8)  Y(co,co')   =  |  (co)  6(co-<o»)   . 

Equation-  (8)  is  sometimes  interpreted  by  saying  that  in  the  case  of  stationary 
random  processes  the  two-dimensional  spectral  density  ¥(&!,«>•)  can  be  localized 
on  the  line  co  =  co'  in  the  (co,co') -plane,  because  of  the  properties  of  the 
delta  function,  (8)  can  just  as  well  be  written  as 

(8)  Y(co,co«)   =   1^^)  6(co-co«)   j 

the  significance  of  this  form  will  appear  below. 


KJ,[2] 
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2.  Locally  stationary  random  processes 

Locally  stationary  random  processes  are  a  natural  generalization 
of  stationary  random  processes,  and  will  be  defined  as  follows  .  The  zero- 
mean  random  process  x(t)  is  said  to  be  locally  stationary  (in  the  wide  sense) 
if  its  covariance  can  be  written  as 


(10) 


r<t,f>  -  <|x(L^A|2>  net-to  , 


where  as  before  R(t-t')  is  a  normalized  correlation  function.  The  factor 
<^  x  ( — * — )    y    represents  change  of  the  average  instantaneous  power  with 
time;  we  shall  assume  that 


I 


|x(r)|2  dT  <  oo 


except  for  the  special  case  <^|x(T)~|^>  ■  const.,  corresponding  to  station- 

t  + 1  • 

ary  processes.  The  symmetrized  variable  — * —   has  been  chosen  to  assure 

that  ]   (t,f)  =  j   (t',t)  as  it  must,  and  because  of  its  suggestive  mean- 
ing as  the  average  of  the  points  t  and  t1.  Heuristically  speaking,  equation 
(10)  should  be  a  good  approximation  if  the  average  instantaneous  noise  power 
changes  slowly  with  respect  to  the  characteristic  width  of  the  correlation 
function  R.   For  simplicity,  we  call  the  functions  appearing  in  the  right 
side  of  (10)  y?    and  L1  respectively,  and  rewrite.  (10)  as 

(id  ru.f)  -  (?  A-^l  it  (*-*■>  • 


(T^) 


In  this  connection,   see  the  definition  of  local  homogeneity  given  by 

■ca. 


A.N.  Kolmogorof f  L  J  . 
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We  now  come  to  the  main  point  of  this  note,  which  is  that  if  a 
random  process  x(t)  is  locally  stationary,  so  that  its  covariance  \      (t,f) 
is  of  the  form  (11),  then  the  same  is  true  of  its  two-dimensional  spectral 
density  T(co,oi')»  To  see  this,  we  first  write  (5)  in  the  form 

!(»,«»)   =  -L^  j   J  expUfe^)  (t-t')-i(^-)  (^')ln(t,t')dtdt»  . 
Then,  introducing  the  new  variables  u  =  — * — ,  v  ■  t-t'  and  using  (11), 


we  find  that 


H2/ 


(12)  Y(»,»')     =     ^il^^-j     V"-"') 

as  asserted,  where 


(13)      v-)    =  & 


exp(-ivw)  fT  (v)  dv  , 


V 


>(z)  =  2^    exp(-iuz)  j^1  (u)  du  , 
and 
(1U)  £  M   =     exp(iuz)  T2(z)  dz  , 


[^  (v)  =     exp(ivw)  Y1(w)  dw  . 


Thus 


the  functions  I,  and  f?  as  well  as  Y?  and  [7*  are  Fourier 


transform  pairs.  It  follows  that  Y.  is  non-negative,  being  the  Fourier  trans- 
form of  the  correlation  function  ["?  ,  and  that  Y~  is  a  correlation  function, 
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being. the  Fourier  transform  of  the  non-negative  function  [?  •  Note  that 
when  Y?  is  the  delta  function,  (12)  reduces  to  the  spectral  density  of  a 
stationary  random  process  (see  equation  (9)). 

Equations  (11),  (12),  (13)  and  (Ik),   together  with  the  stated 
requirements  on  the  functions  [71  ,  |T  ,  Y, ,  and  Y?  constitute  a  natural 
generalization  of  the  usual  Wiener-Khintchine  relations  for  stationary 
random  processes.  Moreover  our  study  of  locally  stationary  random  pro- 
cesses has  "  revealed"  the  structure  of  stationary  random  processes. 
The  delta  function  in  equations  (8)  and  (9)  for  the  two-dimensional 
spectral  density  of  a  stationary  random  process  is  the  Fourier  trans- 
form  of  the  constant  average  power  <^|x|  ^>,  just  as  the  correlation 
function  of  stationary  white  noise  (with  constant  power  spectrum)  is 
the  delta  function. 

3.  Conclusions 

The  concept  of  a  locally  stationary  random  process  introduced 
here  may  prove  of  value  in  technological  applications  of  noise  theory, 
since  local  stationarity  is  a  much  more  reasonable  assumption  to  make 
about  the  kind  of  random  processes  encountered  in  nature  than  is  station- 
arity. (However,  the  usual  heuristics  of  approximating  locally  stationary 
random  processes  in  different  regions  by  "pieces "  of  suitable  stationary 
random  processes  is  undoubtedly  sound.)  Consider,  for  example,  the  pro- 
blem of  the  scattering  of  radio  waves  by  dielectric  noise,  i.e.,  by 
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turbulence-induced  fluctuations  of  refractive  index.  If  turbulent  mixing 
operating  in  gradients  of  refractive  index  is  accepted  as  the  mechanism 
responsible  for  the  dielectric  noise  (as  described  elsewhere  ^-l  *  IAI  >  LPj )  ^ 
then  it  is  clear  that  the  power  level  of  the  dielectric  noise  will  vary 
in  space  and  time  in  proportion  to  the  prevailing  local  gradients  of  re- 
fractive index.  On  the  other  hand,  aside  from  this  variable  power  factor, 
the  local  structure  of  the  turbulence  can  be  expected  to  be  the  same  every- 
where, at  least  if  attention  is  confined  to  a  range  of  eddy  sizes  intermediate 
between  the  macro-eddy  and  micro-eddy  dimensions.  This  situation  is  typical 
of  those  to  which  the  notion  of  locally  stationary  random  processes  seems 
well  suited,  and  indeed  it  is  the  physical  problem  which  suggested  the 
notion  to  the  author. 

The  mathematically  inclined  reader  will  find  some  material  rele- 
vant to  this  note  in  papers  by  Getoor'-  J  and  Loeve'--'. 
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